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Abstract 



(N 

^ ' We establish existence and uniqueness of generalized solutions to the initial-boundary value 

problem corresponding to an Euler-Bernoulli beam model from mechanics. The governing par- 
tial differential equation is of order four and involves discontinuous, and even distributional, 
coefficients and right-hand side. The general problem is solved by application of functional 
analytic techniques to obtain estimates for the solutions to regularized problems. Finally, we 
04 ' prove coherence properties and provide a regularity analysis of the generalized solution. 
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- - ■ 1 Introduction and basic estimates 

We consider the question of existence and uniqueness of a generalized solution to the initial- 
boundary value problem 

u|t=o = fi,dtu\t=o = /2, (1) 
mU=o = = 0, dxu\x=() = dxu\x=i = 0, 

where Q is a differential operator of the form 

Qu := dl{c{x)dl)u + b{x, t)dlu, 

and b,c,g,fi and /2 are generalized functions. The precise structure of the above problem is 
motivated by a model from mechanics describing the displacement of a beam under axial and 
transversal forces, which we briefly present in Subsection 1.1. 
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In Subsection 1.2 we review basics from Colombeau's nonlinear theory of generalized functions, 
which are the framework of our main results in Section 3. 

In Subsection 1.3 we briefly sketch the well-known functional analytic tools from the variational 
approach and give a refined version of energy estimates, which lie at the heart of the proofs in 
Sections 2 and 3. 

Section 2 is devoted to existence and uniqueness of solutions in the Sobolev space setting for 
(slightly) regularized coefficients b and c. 

Finally, Section 3 establishes existence and uniqueness of a solution to the original problem 
in the space of Colombeau functions based on families of weak solutions to the corresponding 
regularized problem. Furthermore, we prove coherence of the generalized solution with classical 
and weak solutions in case of smooth coefficients but also analyze generalized regularity properties 
in the non-smooth case. 

1.1 The Euler-Bernoulli model 

This subsection serves as a motivation for the analysis of problems of the type H]). Consider the 
Euler-Bernoulli beam under a distributed transversal force g and an axial force P (cf. [2]). As has 
been initiated in [16j we pursue to investigate the case when the beam consists of two parts with 
different physical properties and cross sections. Let u denote the displacement depending on the 
spatial variable x and time t. Accordingly the differential equation of the transversal motion reads 

where: 

• A denotes the bending stiffness and is given by A{x) — EIi + H{x — xo)El2 — EIi. Here, 
the constant E is the modulus of elasticity and Ji, I2 {Ii ^ I2) are the moments of inertia 
that correspond to the two parts of the beam and H is the Heaviside jump function. 

• R denotes the line density (i.e., mass per length) of the material and is of the form R{x) — 
R„ + H{x-xo){Ri-R2). 

• P represents the axial force, given by a time dependent function or distribution. For example, 
it may be of the form P = Pq + PiS{t - to) where Pq > 0, -Pi > and to > 0, or P = 
Po -1- Pi sm{ujt) with Po > 0, Pi > and oj > 0. 

' g is the force term and can be of the form g{x, t) — Fq5{x — xi), where xi € (0, 1). Further- 
more, for applications also a specific dependence of g on u is of interest, as, for example, in 
the case of the so-called Winkler foundation with g = —cu, c > constant. 

• As mentioned above, u = u{x,t) denotes the displacement. Thus, its second derivatives 
and represent acceleration and linearized curvatur({3, respectively. 

In addition to ([2]) we consider initial conditions: 

w(x,0) = /i(x), dtuix,0) ^ f2ix), (3) 

where fi and /2 are the initial displacement and the initial velocity. We consider the beam to be 
fixed at both ends and also supply boundary conditions of the form 

u(0,t) = w(l,t) = 0, d^u{0,t) = dM^,t) =(^- (4) 

By a change of variables t t via t{T) — y/ R{x)t we transform the problem ((^E]) into the 
standard form given in ^ . The function c in ([T]) equals A and therefore is of Heaviside type and 

^without linearization it would read k = j-. 
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the function b is then given by b{x, t) ~ P{R{x)t) and its regularity properties depend on the 
assumptions on P and R. 

We show in Section 2 that standard functional analytic techniques reach as far as the following: 
boundedness of c and 6 together with sufficient (spatial Sobolev) regularity of the initial values 
fi, f2 as well as of g ensure existence of a unique solution u G L^{0, T; H^{{0, 1))) to ([1]). However, 
the prominent case P = Pq + PiS{t — to) is clearly not covered by such a result. Moreover, since 
both c and d^u are distributions their product is not always defined. In order to allow for these 
stronger singularities one needs to go beyond distributional solutions. 

1.2 Nonlinear theory of generalized functions 

We will set up and solve problem ([1]) in the framework of algebras of generalized functions on 
the domain Xt = (0, 1) x (0,T) (with T > 0). To be more specific we employ a variant of the 
Colombeau spaces based on L^-norm estimates as introduced in [3] (and applied in [T^). Needless 
to say that the general theory is set-up on arbitrary open subsets of R"* or even on smooth 
manifolds, but we use this opportunity to introduce notation appropriate for our specific model 
(cf. |5l [51 El US]). The basic objects defining our generalized functions are regularizing famihes 
('^e)ee(o,i] of smooth functions Ug G H°°{Xt), the space of smooth functions on Xt all of whose 
derivatives belong to L^. To simplify notation we will henceforth write (ite)e to mean (we)6e(o,i]- 
We single out the following subalgebras: 

Moderate families, denoted by £ m ,h°° (Xt) ■• are defined by the property 

Va e NS, 3p > : Wu^Wl^^x^) ^ 0(e-f ), as e 0. 

Null families, denoted J^h'^(Xt)^ are the families in £m,h^{Xt) satisfying 

Vg > : ||ii,|U2(;,^) = 0(e«) as£->0. 

Hence moderate families satisfy estimates with at most polynomial divergence as e — > 0, 
together with all derivatives, while null families vanish very rapidly as e — s- 0. For the latter, one 
can show that, equivalently, all derivatives satisfy estimates of the same kind (cf. [9j Proposition 
3.4(ii)]). The null families form a differential ideal in the collection of moderate families. 
The Colombeau algebra is the factor algebra 

Qh°°{Xt) ^ ^MM'"{Xt)/-^H°°{Xt)- 

We will occasionally use the notation [(ug)^] for the equivalence class in Gh=°{Xt) with repre- 
sentative (ue)e. The algebra Gh°° {(q.i)) on the interval (0, 1) is defined in the same way and its 
elements can be considered being members of Qh^(Xt) as well. 

Remark 1.1. (i) For any (ue)e in £m,h°^{Xt) have smoothness up to the boundary, i.e. € 
C°°([0, 1] X [0,T]) and therefore the restriction u\t=o of a generalized function u G Gh°^{Xt) 
to t = is well-defined by u^{-,0) £ £m.h°°((q,i))- (As emphasized in |31 Remark 2.2] this 
follows from general Sobolev space properties, since Xt clearly is a Lipschitz domain; cf. 

m.) 

(ii) If in addition to w G Qh°°((o,i)) we have for some representative (we)e of u that G Hq{{0, 1)), 
then Vg{0) — Ue(l) — and dxVs{0) — dxVe{l) = (cf. Remark |2. II below) . In particular, 

v{0) = v{l) = and dxv{0) = d^vil) = 

holds in the sense of generalized numbers. 

(iii) Note that i'^-estimates for parametrized families S H°°{Xt) always yield i°°-estimates of 
the same qualitative behavior with respect to e (since H°°{Xt) C C°°{Xt) C W°°'°°{Xt))- 
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As described in [7 the space H~°°{M.'^), i.e. distributions of finite order, are embedded (as a 
Unear space) into Qfj^^^d-^ by convolution regularization. Furthermore, by this embedding H°°(R'^) 
becomes a subalgebra of GH°°(m)- 

Some generalized functions display so-called macroscopic or distribution aspects: We say that 
u = [{ue)e] G Qh°° is associated with the distribution w € V , denoted hy u ^ w, if for some (hence 
any) representative {ue)e of u we have Ue ^ w in V , as £ ^ 0. 

Intrinsic regularity theory for Colombeau generalized functions has been started by defining 
(in T^) the subalgebra Q°° , which satisfies the crucial compatibility property Q°° 2?' = C°° . In 
this sense, it plays the same role for Q as C°° does for V . Similarly 13J introduces the subalgebra 
Q'^ao of regular elements in Qh°° by the following condition: u G Qh-^^ belongs to C/^oc. if and only 
if 

3p > 0,Va e NJJ : Wu^\\l-- = 0(£~p), as e ^ 0. (5) 
Note that p can be chosen uniformly with respect to a. 

1.3 Precise constants in energy estimates from variational methods 

Let y, H be two complex, separable Hilbert spaces, where V is densely embedded into H . Denote 
the norms in V by | • | and in by || • ||. Thus, if V' denotes the antidual of V , then V C H C V 
forms a Gelfand triple. 

Let a(t, .,.), ao(t, .,.,), and ai(i, .,.) (t € [0,T]) be (parametrized) families of continuous 
sesquilinear forms on V with 

a{t,u,v) = ao{t,u,v) + ai{t,u,v) Vu, u G (6) 

such that the "principal part" oq and the remainder oi satisfy the following conditions: 

(i) for all u,v £ V: < i— s- ao(t, u, v) is continuously differentiable [0, T] C, 

(ii) flo is Hermitian, i.e. ao(t, u, v) — ao{t, v, u) for all u,v V, 

(iii) there exist real constants A and a > such that 

ao{t,u,u)> a\u\^ - X\\u\\^ e V, Vt e [0,T], (7) 

(iv) for all u,v €i V: 1 1-^ ai{t, u, v) is continuous in [0, T] — s- C, 

(v) there is Ci > such that for all t e [0,T] and u,v eV: \ai{t,u,v)\ < Ci\u\ \\v\\. 

Note that writing a'Q{t,u,v) = ^(ao(t,u, u)) condition (i) implies the following: there exist 
nonnegative constants C and Cq such that for all t e [0, T] and u,v G V 

\ao{t,u,v)\ <C\u\\v\ and \a'Q{t,u,v)\ < Cq\u\\v\. (8) 

As shown in 8, Chapter XVIII] the above conditions guarantee unique solvability of the ab- 
stract variational problem as described by the following Theorem. 

Theorem 1.2. Let a{t,.,.) (t e [0, TJJ satisfy conditions (i)-(v). Let uq E V , ui E H. and 
f E i^((0, T), H). Then there exists a unique u E i^((0, T), V) satisfying the regularity conditions 



u'=^EL\iO,T),V) and = J G £^((0, T), O (9) 

and solving the abstract initial value problem 

(u"{t),v) +a{t,u{t),v) (f{t),v) yvEV,ytE{0,T) (10) 

u(0)=uo, u'(0) = ui. (11) 



(Note that (0) implies that u E C{[0,T],V) and u' E C([0,T],y). Hence it makes sense to 
evaluate m(0) E V and u'(0) E V' and claims that these equal uq and ui, respectively.) 
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The prove of this theorem in [51 Chapter XVIII] (see also [T71 Chapter III, Section 8]) proceeds 
along the following lines: first, one shows that u satisfies a priori (energy) estimates which imme- 
diately imply uniqueness of the solution and then deduces existence of a solution by a Galerkin 
approximation. To identify the precise dependence of all constants in the a priori estimates we 
adapt the corresponding part of the proof to our case (which happens to be simpler than in [8] , 
whereas |17j does not cover our situation). 

Proposition 1.3. Let u be a solution to the abstract variational problem [ WTT]) . then for all 
t€[0,T] 

\u{t)\' + \\u'{t)r < {Dt \uo\' + Ibiif + r \\f{r)fdT) ■ eMtFr) (12) 



where the constants Dt and Ft are given in terms of the constants C,Co,Ci,a, X occurring in 
conditions (iii),(v), and ([8]) explicitly by 

= (C + A(l +r))/min(a,l) and Ft = max{Co + Ci, Ci + T + 2}/ min(a, 1). 

Proof. We put v = u'{t) in (fTU]) and obtain 

{u"{t),u'{t))+a{t,u{t).n'{t)) = {f{t),u'{t)). (13) 

By sesquilinearity and (ii) we have 

|h'Wf - |((u'(t),H'(t)» = {u"{t),u'it)) + («'(<), ^^"(^)) = 2Rei(u"it),u'm 



and a{t, u, v) = ao(i, v, u) + ai(t, u, v). Thus taking real parts in (fT3)) we obtain 

^\\u'{t)f + ao(i, u'{t)Mt)) + ao(i, u{t), u'{t)) 
dt 



+ai{t,u{t),u'{t))+ai{t,u{t),u{t)) = 2Re(/(t), u'(t)). (14) 
For any v e H^{{0,T),V) 

^ {ao{t, uit),v{t))) - a'oit, u{t),v{t)) + ao(i, u'it),v{t)) + aoH, uit), v'it)), 
hence implies 

^{\\u'{t)\\^ + ao{t,uit),u{t))j = -2Reai{t,u{t),u'{t))+a'„it,u{t),uit)) + 2Re{fit),u'it)) 

(15) 

Integration (fT5)) we obtain 

LHS{t) ■.=\\u'{t)f + ao{t, u{t), u{t)) 

= ||m'(0)|P + ao(0,u(0),M(0)) - 2 / Reai(T, u(r), u'(r))(iT 



+ / ao(r,u(r),u(r))dr-|-2 / Re {f{T),u' {T))dT =: RHS{t). (16) 



Further, using the conditions (i) and (v), initial conditions ([TT|) and the inequality 2ab < + b^ 
we obtain that the right-hand side of (flB]) can be estimated in the following way: 



\RHSit)\<\\uif + C\uof + 2CiJ \u{T)\\\u\T)\\dr + Co I \uiT)fdT 

+ f\\f{T)fdT+ f\\u'{T)fdT 



<C|uoP + hi|P + (Co + Ci) / \u{r)\'dT+ I |l/(r)f dr + (l + Ci) / \\u'{T)fdT. 

(17) 
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Condition ([7]) implies that the left-hand side of has the lower bound 

LHS{t)>a\u{t)\-X\\u{t)r + \\u'{t)r, 

and therefore 

a\\u{t)f + \u'{t)f < \RHSit)\ + XHt)\\ (18) 
We claim that for all te [0, T]: 

\u{t)\^<{l + minor + f\u'{T)\dT). (19) 



Indeed, since u{t) = uq + /p u'{T)dT and ||uo|| < |mo| we have 



\\u{t)f ^ {u{t),u{t)) =\\uor + 2^0 [ {uo,u'{T)dT) + { [ u'{T)dT, [ u' {s)ds) 

Jo Jo Jo 

<||U0||"+ j'lluafdT + j \\u'{T)\\'dT + j j \{u'{T),u'{,))\d,dT 

Jo Jo Jo Jo ^ 

< ||uo||'+i||uoir+ / \\u{T)fdT + t f \\u'{T)\\dT 

Jo Jo 

<(l + i)(lko|P+ /V'(r)|pdr)<(l + t)(KP+ f\\u'{T)rdT). 
Jo Jo 

Let [3 := min{Q;, 1}. Combining ([TT]), and ^5)1 we arrive at 

/3(|u(t)p + h'(t)lp) < (C + A(l + i))|"oP + ll^^ilP + r ||/(r)|P dr+ 

Jo 

(Co + Ci) / |w(r)|2dr + (2 + t + Ci) / \\u'{T)\\dr 
Jo Jo 

T 

2 I IL, l|2 I / II j/„M|2, 



< (c + A(i + ^))|^ior + hiir + / ii/wir^T 

+ max{Co + Ci,Ci+r + 2} / {\u{t)\' + \\u' {T)f)dT. 

Jo 

Dividing by /? > yields 

m<DT\uof + \\uif+[ \\f{T)fdT + FT [ <^{T)dT, 

Jo Jo 

where $(i) = |u(i)|2 + |ltt'(t)||2 (^^^^j ^ (C + A(l + r))//3 and Ft = max{Co + Ci, Ci +T + 2}//3 
as in the statement of the Proposition). Gronwall's lemma now implies that 

m < {Dt \uo\^ + lluilP + r ||/(T)f dr) • eMtFr). 

Jo 

■ 

2 Weak solutions for L°°-coefRcients 

Let H -^^((0, 1)) with the standard scalar product {u, v) = u{x)v{x)dx and L^-norm denoted 
by II '11. Let V be the Sobolev space Hq{{0, 1)), which is the completion of the space C^((0, 1)) 



(compactly supported smooth functions) with respect to the norm ||u||2 = (X]fe=o H'"'''^-' IP)'^^^ (and 
inner product {u,v) i-^ X]fc=o ''^''^''))- Then V' = H~^{{0, 1)), which consists of distributional 
derivatives up to second order of functions in L^{{0, 1)), and V ^ H ^ V' forms a Gelfand triple. 

Let c e i°°((0,l)) be real-valued and b e C([0, T], L°°((0, 1))). For t e [0,T] we define the 
sesquilinear forms a{t, •, •), ao(t, •, •), and ai{t,-,-) on x by 

ao{t,u,v) ^ {cdlu,dlv), ai{t,u,v) ^ {h{t)dlu,v) {u,v eV) (20) 

and 

ait, u, v) = ao(t, u, v) + ai{t, u, v). (21) 

Clearly, oq is Hermitian (since c is real). 

Note that the L°°-properties of c and b are necessary in order to have the above sesquilinear 
forms defined and continuous on all of V. In the application to our model (jT]) the functions c and 
b play the role of the coefficients (with the abuse of notation in identifying b{t){x) with b{x,t)) 
and hence the required L°°-condition is satisfied for c but typically not for b, which we want to 
allow to be a measure or more general distribution. This poses the main difficulty in the analysis 
of the model. Our strategy will be to regularize the problem by smoothing of the coefficients 
and data and then analyze the corresponding family of solutions and show that they constitute a 
generalized solution. 

As a preparatory result, we need precise statements and estimates for the solutions to the 
regularized problems. In fact, for the basic result we will have to assume considerably less regularity 
than actually needed in the theory applied later on. The advantage is that it illustrates how far 
we could go by classical functional analytic methods concerning lowest possible regularity of the 
coefficients. At this stage we make the hypotheses 

cei°°(0,l), 5(E C([0,r],L°°) (22) 

and (motivated by the specification of the bending stiffness in Subsection 1.1) that there exist 
Ci > Co > such that 

< Co < c{x) < Ci (for almost every x). (23) 

Remark 2.1. (About boundary conditions) Note that any solution of the partial differential 
equation and the initial conditions in ([1]) with u{.,t) E Hq{{Q,1)) for all t e [0,T] automatically 
implements the boundary conditions of ([1]) in the following sense: -ffo((0, 1)) is continuously 
embedded in the subspace of C^([0, 1]) of functions with vanishing values and derivatives at the 
boundary ([^D, Corollary 6.2]); hence we must have u{0,t) — u{l,t) = as well as dxu{0,t) — 

Theorem 2.2. Let b,c be as in i f^) and I12S\) and the sesquilinear forms a(t, .,.) (fort G [0,T]) 
be defined by ^ and If fi G H^{{0,1)), /2 G i^((0,l)), and g € L^{{0,T), L^{0,1)), then 

there exists a unique u S L^((0, T), TJq ((0, 1))) satisfying 

and solving the initial value problem 

{u"it),v)+a{t,u{t),v)^{git),v), Vi; G ij2((0, 1)), t G (0, T), (25) 
u{0) = fi, m'(0)=/2. (26) 

(Regarding the precise meaning of the initial conditions 1126]) we refer to the corresponding remark 
in Theorem ) 

Proof. We show that the sesquilinear form a satisfies all hypotheses of Theorem 11.21 First, 
for all M, u G i7o((0, 1)) we clearly have that oq is Hermitian (since c is real) and continuously 
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differentiable (by independence of t) and that ai is continuous with respect to t and satisfy the 
estimates 

\ao{t,u,v)\ < ||c|U=o((o,i))p2^||PM <C\\uUvh 

\ai{t,U,v)\ < |j6||Loo((o,l)x(0,T))||^^w||||«|| < Ci ||u||2||w||2, 

where C = |lc||i=o((o.i)) and Ci = ||6||l~((o,i)x(o,t))- Hence conditions (i), (ii), (iv), and (v) are 
met. 

To prove property (iii) we appeal to EhrUng's lemma ([501 Theorem 7.4] or [4, Chapter 2, 
Proposition 2.11]) which implies that for each S > there exists a real constant Cs such that for 
ahue ff2((0,l)) 

\\u\\l<S\\u\\l + Cs\\u\\^ (27) 
Choosing S ~ 1/2 and employing ([23]) we obtain for all u e iJ^((0, 1)) 

ao{t, u, u) — {cd'^u, d^u) > co||9^w||^ = co||w||2 — cq ||w||i 

> coll^ll^ - coi^WnWl + Ci/2h|P) > y - coC./^M^ 

= a\\u\\l-\\\ur, 

where a = cq/2 and A — C1/2C0. Thus we have shown also (iii). ■ 

Remark 2.3. (i) For later reference we give the precise dependence of all constants appearing 
in the energy estimate (|12p 

\\uml + \\u'it)f < [Dt holl^ + ll^^ilP + r Mr)fdT) ■ expitFr) 

Jo 

on the coefficients b and c. We recall that 

Dt = (C + A(l + r))/min(a,l) and = max{Co + Ci, Ci + T + 2}/ min(a, 1), 

where we now have 

C = I|c||l°o((o^i)), Co = 0, Ci = ||&||l=°((o4)x(o,t)) 
a = y, A = Ci/2 Co- 
in) According to jU) the solution u belongs to C^{[0,T], H-^{{0, 1))) ^ V'{{0, 1) x (0,r)) and 
in case of smooth coefficients b and c is a distributional solution to the partial differential 
equation 

d^u + dl{cdl)u + bdlu^g inl?'((0,l) x (0,r)). 

Moreover, if uq, ui, and g are smooth, then m is a classical smooth solution to the partial 
differential equation. (We prove this in course of the proof of Theorem l3.1l below: cf. Remark 

|3j;iii).) 

(iii) In view of the statement in (ii) it might be interesting to note that the partial differential 
operator in our model does not fall into any of the standard types of linear differential 
operators of distribution theory. To make this precise it suffices to consider the case of 
constant coefficients c > and & G M. Thus we have the operator 

P = d^ + cd1^ + bdl 

with symbol p{^,t) = — + c^^ — b^^ and principal symbol ^'2(^1''') — c^^. Then the 
following holds (for clarity, with slight logical redundance in the statements): 
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(a) Clearly, P is not elliptic (T, p. 176]), since Char(P) = {{£,,t) G \ {0} : P2(C,t) = 
0} = {(0, r) : T ^ 0} is not empty. 

(b) P is not parabolic (with respect to forward time, i.e. with respect to iJ = f > 0; 7, pp. 
202-203, Corollary 2]), since with ^ > sufficiently large such that t(0 := c^-^-fe^^ > 
we have |^| + |r(^)| ^ og {£, ^ oo), while p{t{() ■ (0, 1) + 0)) = 0. (Note that, in 
contrast, dt + — is parabolic.) 

(c) P is not hyperbohc (with respect to forward time; J7, pp. 190-193, Prop. 16]), since 
(0,1) e Char(P). 

(d) Moreover, P is not hypoelliptic (thus cannot be parabolic or elliptic), since the equiv- 
alent condition [3 (4.11), p. 233] can be shown to be false (cf. [3 Theorem 1 on p. 233 
below]). With k = (0, 1) (i.e. the partial r-derivative) and ^ > sufficiently large to 
ensure r(^) := c^"' — 6^^ > 1 we obtain 19rp(C,T(^))] = 2t(^) — > cx) (as ^ — > oo), 
whereas p(^, t(^)) = for large ^ > 0. 

(e) Of course, P is also not of Schrodinger type (in the sense of [S] pp. 620]). However, 
if = one can write P = {dt ~ i^&^j^ ■ (pt + i\/cd^), which is a product of two 
Schrodinger operators. 

3 Colombeau generalized solutions 

We now establish existence and uniqueness of a generalized solution to problem ([T]), where the 
coefficients 6, c as well as the data /i, /2 and g are Colombeau generalized functions. This means 
that we find a unique solution u G Qh'^{Xt) ^hc partial differential equation 

d^u + Q{t,x,d^)u = g on Xt = (0,1) x (0,T), 

where g G Gh=°{Xt) Q denotes the differential operator on Gh°°{Xt) which acts on represen- 
tatives by 

(Ue)e ^ {dUc{x)dl{Ue))e+b{x,t)dl{u,)), =: Q^Ue, 

with initial conditions 

u\t=Q = /i G 5//°° ((0,1)), dtu\t=o = /2 e Gh°^({o,i))- 

Recall that the initial conditions are to be understood in the sense of Remark ILlT i') . 
Furthermore, the boundary conditions read 

u\j:=0 = u\x=i = 0, dxU\x=0 = dxU\x=l = 0. 

Thanks to Remark lLlT ii) these are automatically satisfied if we can show that for all t G [0, T] we 
have u^{t) G Hq{(Q, 1)) for some representative {ue)e of u. 

As in the classical specification with we have to impose a compatibility condition concering 
initial and boundary values 

/i(0) = /i(l) = (28) 

(which is an equation in generalized numbers). We observe that for any /i G Gh°^{{o,i)) condition 
(PS)) implies that there exists a representative {fi.e)e of /i such that /i^^ G Hq{{0,1)) for all 
e g]0, 1[. (If {fe)e is an arbitrary representative, put :— /i,e(0), :— /i,e(l) and consider the 
new representative fi,six) := fi,e{x) — — (mg — Ue) x instead.) 

3.1 Existence and uniqueness 

As in Section 2 we impose a condition on c, which is motivated by the intended properties of 
the bending stiffness: There exist real constants ci > cq > such that c G Gh°° (o,i)) possesses a 
representative {ce)e satisfying 

< CO < ce{x) < ci VxG (0,l),VeG]0, 1]. (29) 
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(Note that then any other representative {ce)e of c satisfies ^ — Ce{x) < ci + 1 for all x G (0, 1) 
and < e < So with some Sq g]0, 1].) A weaker condition on c is briefly discussed in Remark 
inUi) below. 

As a technical assumption we will also require that b is of L°°-log-type (similar to 18 ), which 
means that we have for some (hence any) representative (b^)^ of b there exist N ^ N and eo e]0, 1] 
such that 

||6,|Uoo(x,)<7V.log(i) (0<e<eo). (30) 

As observed in [18, Proposition 1.5] log-type regularizations of distributions are easily obtained 
via convolution with lograithmically scaled mollifiers. 

Theorem 3.1. Let b G Gh'^(Xt) '^f -log-type and c G GH-^io.i)) satisfy (j29p . For any 

fi € Gh=°{(o,i)) satisfying (|28| . /2 G Gh°°({o,i)), o,nd g G Gh°^{Xt) there is a unique solution 
u G Gh'^(Xt) to the initial-boundary value problem (Qp. 

Proof. Existence: Let {b^)e represent b and (ce)^ be a representative of c satisfying (jSO]) . 
Denote by (/ie)e, (/2e)e, and {ge)e representatives of /i,/2, and g, respectively. In addition, we 
may assume that fi,^ G Hq{{0, 1)) for all e g]0, 1[ (see the discussion following 

Let e G]0, 1] be arbitrary. By Theorem 12.21 and Remark I2.3r ii) we obtain a unique function 
Ue G H\{0, T), H^iiO, 1))) n H^{{0, T), i/-2((0, 1))) which solves 

PeU^ := dfU^ + Q^{t, X, dx)ue = g^ on Xt, (31) 

Ue\t=0 = /le, dtUe\t=0 = f2e- 

In particular, we have G C^i[0,T], H-^{{0,1))) n C{[0,T], H^{{0,1))). 

Moreover, from Proposition 11.31 and Remark l2.3r i) we deduce the energy estimate 



ml + IKWIP < {D^T WfuWl + ll/2,e|P + r heirW dr) ■ exp(t F|) 



(32) 



where with some N we have for small e > 



Dt = (llcelk- + A(l + r))/min(a, 1) = 0{\\ce\\L^) = 0(1) (33) 
F| = max{||&,||ioo, \\b,\\L^ +T + 2}/min(a, 1) - OQMl^) ^ 0(log(e-^)), (34) 

since a — co/2 and A = C1/2C0 are independent of e. Therefore moderateness of the initial 
conditions fie, f2e and of in ([32|) yields that there exists M such that for small e > 

Ikelli.fx.) + + WdlueWl^^^^^ + = 0(6-'') {e ^ 0). (35) 

We will proceed to show that the family (we)e belongs to £h°°{Xt)- Then by construction its 
class u in Gh°°(Xt) defines a solution to the initial value problem. It remains to prove the following 
properties: 

1. ) For all e G ]0, 1] the function is smooth, i.e. G C°"{Xt)- 

2. ) Moderateness, i.e. for all fc G N there is some A/ G N such that for small e > 



Note that <^ already yields for (?, k) G {(0, 0), (1, 0), (0, 1), (0, 2)}. 

Step 1: Differentiating (considered as an equation in I?'((0, 1) x (0,T))) with respect to t we 
obtain 

P,{dtu,) = dtge - dtb{x,t)dlu, =: ~g, G H\{{),T),L\{Q,1))), 



10 



since dtge e H°°{Xt), dtb{x,t) G H°°{Xt) C W^'^{Xt) and dlu^ G H^{{0,T),L^{{Q,l))). 
Furthermore, since depends smoothly on t as a differential operator in x and Ue(0) — fi^^ G 
iJ°°((0, 1)) we have 

(9iU,)(.,0) = /2,e=:/i,eei^°°((0,l)), 

{dt{dtuM;0) = ffe(.,0) - Qe(«e(.,0)) = ffe(.,0) - Q,/i,, h,, G i/°°((0, 1)). 

Hence (9fite satisfies an initial value problem for the partial differential operator as in pip 
with initial data /i.^, /2,e and right-hand side instead. However, this time we have to use 
V = H^{{0, 1)) (replacing H^{{0, 1))) and H = L^{{0, 1)) in the abstract setting, which stiU can 
serve to define a Gelfand triple V ^ H ^ V' (cf. [2Ql Theorem 17.4(b)]) and thus allows for 
application of Theorem 11.21 and the energy estimate (fT2|) (with precisely the same constants) . 

Therefore we obtain a* G H'^{[0,T],H^{{0,1))), i.e. G H^{{0,T), H^{{0,1))) and from the 
variants of (|32p (with exactly the same constants Dj. and Fj.) and ((35| with 5fMe in place of 
that for some M we have 

+ \\d..dtu,\\l.^^^^ + \\dldtu,\\l.^^^^ + = 0(£-*') (e - O). (36) 



Thus we have proved (TJ^fc ) with (/, k) = (2, 0), (1, 1), (1, 2) in addition to those obtained from (I35p 
directly. 

Differentiating pip a second time with respect to t we obtain again an initial value problem 
for the same operator with solution dfu^: this time, the right-hand side reads 

I, := dfg, ~ d1h{x,t)dlu, - dtb{x,t)dldtu, 
and the initial values are /i^e := {dfU^){.,0) — f2,e and 

he ■■= {dtid^tuM;^) = at.9e(.,0) - Qeike) " dM.,0)dlh,,. 

Note that the corresponding energy estimate involves again precisely the same constants Dfp and 
as in p2p and we only have to insert the norms of the new initial values and right-hand side 



accordingly. Thus we obtain G ((0, T), ((0, 1))) and L -norm estimates proving (Ti^k 
also for the new pairs (/, k) = (3, 0), (2, 1), (2, 2) in addition. 

Upon deriving a succession of similar initial value problems for dfu^, dfu^ etc. we arrive at 

u,eH°-iiO,T),H'iiO,m (37) 

and 



^Ilk\ holds for alU G N and fc = 0, 1, 2. (38) 
Using this information in the original equation pip gives 

Q,u, ^g,- d^u, G i/°^((0, T), i/2((o, 1)). 

Since is an elliptic operator of order 4 (in the a;- variable), p7p and elliptic regularity implies 
Me G i?°°((0, T), 7J^((0, 1))). Using this improved regularity property in turn we obtain by a 
successive application of elliptic regularity 

ue e h^{Xt) c c°°(x^). 

Thus, requirement 1.) above is proved. 

It remains to show the moderateness property 2.). Due to (I38p the statement is already verified 
for derivatives of arbitrary orders with respect to time and up to order 2 with respect to x. 
Step 2: From the results of Step 1 and equation (|3T|) we deduce that 

he := dl{ce dlue) ^ ge~be dlue - dfUe 
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satisfies for alH G N with some Ni an estimate 

\\dih,h2(x^j ^ 0{e-^') (e->0). (39) 
Integrating the equation ~ d^{ced^Ue) — twice with respect to x gives 

l-x ry 

e^{t) + d^{t)x^Ce{x)dlu^{x,t)- I / h^{z,t)dzdy. (40) 



Evaluating this equation at a: = and a; = 1 we obtain ei,{t) — c^{Q) d'^u^{0,t) and d^{t) — 
Ce{l)dlue{l,t) ~ h^{z,t)dzdy. Thus §^ and ^ show that (e^)^ and (4)^ belong to 

£m,H'^ {{0,1}) ■ 

Equation ((40|) yields the following formula for all I, k £ N, k > 2: 

d\dlu,{x,t) = dl-^ ( d'Mt)+ d[d,^t)- +J- r r d\hAz.t)dzdy) . (41) 

\ Ce[X) Ce[X) Jq Jq / 

Using the fact that is moderate (by boundedness away from 0) we immediately obtain from 
([39)1 and (|4T|l the moderateness estimates for ||9j9^Uj||i2(x^) with fc = 3,4 as well (l arbitrary), 



thus {Tik\ holds for alU e N and fc = 0, . 



Applying now (Tj.fe ) with fc up to 4 in (appropriate derivatives of) the defining equation for 



above we obtain the following improvement of 

V/eN,0<p<2, 3iV,p: WdPdlheWmx^) ^ 0{e-^'-) (e ^ 0). 



Using this in turn in (|4ip implies (T;_fe I for k up to order 6. Successively proceeding in this way, 
we conclude that {Ti^k I holds for all l,k gN. Thus property 2.) is shown and therefore u — [{u^)^] 
defines a solution u € Qh°°(Xt) initial value problem. 

By construction we have u^{t) e Hq{{0, 1)) for all t G [0, T], hence u(0,t) = u{l,t) — and 
dxu{0,t) = dxu{l,t) = and u satisfies the boundary conditions too. 

Uniqueness: Assume that u = [{ue)e] satisfies initial-boundary value problem with zero 
initial values and right-hand side. Then we have for all g > 

= 0(£«), = 0(£'), heWmx^) = O(e^) {s -> 0) 

and the energy estimate (|32|) together with p3ll34p imply for all q > an estimate 
Therefore (Me)e £ J^h^{Xt) which proves that u = 0. ■ 

Remark 3.2. (i) The above proof is valid without change if we replace condition on c e 
Sh°°{(o.i)) by the considerably weaker condition that c is strtictly positive, i.e. there exists 
g > and eo > such that 

c, > e« (0 < £ < £o) (EH) 

(and no extra boundedness from above [other than the standard moderateness]), since this 
clearly implies that 1/c is a generalized function (cf. also jTIj Theorem 1.2.5]). Therefore, 
the existence and uniqueness result is extended to cases with stiffness that may be infinitely 
close to and need not be bounded from above. 

(ii) In Step 2 of the proof we took advantage of having only one spatial dimension. However, 
for a variant in higher space dimensions basic conclusions similar to those we drew from 
formula ((4T|) could be reached by pseudo-differential techniques as in [10] with the help of a 
Colombeau generalized parametrix of the operator v i-^ dxicd^v). 

(iii) We point out that along the way in the proof of Theorem 13. II (namely in Step 1) we have 
proved the following regularity result for the weak solution: 

If in addition to the hypotheses of Theorem 12.21 the coefficients b and c as well the data 
/i, /2, and g are smooth, then the unique solution u is also smooth, more precisely Ue G 
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3.2 Regularity and coherence with smooth and weak solutions 

As a first simple comparison results displaying an enjoyable consistent relationship between clas- 
sical smooth solutions, weak solutions, and generalized solutions we consider the case of C°° 
coefficients. 

Proposition 3.3. Assume in addition to the hypotheses of Theorem \3.1\ that the coefficients b 
and c are (Colombeau classes of) C°° functions. 

(i) If the data fi, /2, and g are (the Colombeau classes of) smooth functions as well, then the 
unique generalized solution u G Gh^{Xt) the initital-boundary value problem ^ equals 
the (embedding of the) unique smooth solution in H°°[Xt)- 

(a) Let fi = [(/i,e)e], /2 — [(/2,e)e]; o,nd g = [{ge)e] 0''^^ assumc that we have the following 
convergence properties of regularized data as e —* 0: fi^e —* fi 7Jg((0, 1)), f2.e /2 in 
L^((0, 1)), and g^ ~* g in L'^{Xt)- Then the unique generalized solution u G Gh°°{Xt) 
associated with the unique weak solution w G i^((0, T), Hq ((0, 1))) according to Theorem \2.S\ 
with initial values fi, f2, and right-hand side g. 

Proof, (i) Put := V for all e, where v denotes the unique classical smooth solution (smooth- 
ness of V is due to Remark 13.2^ 111)). From uniqueness of the Colombeau solution u we must have 

u = [m- 

(ii) This assertion follows immediately from the abstract theorem [SJ Chapter XVIII, Section 5, 
Subsection 4.2, Theorem 2] on continuity of the solution with respect to the initial values and the 
right-hand side. ■ 

Finally, we investigate intrinsic regularity of the Colombeau generalized solution u in the sense 
of i.e. in terms of uniform e-asymptotics for all derivatives. In the context of regularity 

theory for (pseudo-) differential equations a key property required of the coefficients (or symbols) 
has been shown to involve the notion of so-called slow scale nets ( [lOl [Ml [15] ) : A net (r£)5g]Q of 
complex numbers is said to be of slow scale, if 

VpGN: |r^|P = O(e-i) (e -> 0). 

Note that any product of finitely many slow scale nets is of slow scale. A generalized function 
V = [(we)^] G Gh'=° is said to be of slow scale in all derivatives, if the net (||9"we||i2)e is of slow 
scale for all a. 

Assuming slow scale conditions on the coefficients we show that regularity of the initial values 
and the right-hand side is preserved in the solution. This is a generalization of the regularity result 
for weak solutions described in Remark l3.2r iii). 

Theorem 3.4. Let b G Qh=°(Xt) '^'^'^ ^ G Gh°^(o,i)) of slow scale in all derivatives. Assume 
further that c satisfies (|29p and that b is of slow scale L°° -log-type, i.e. for some (hence any) 
representative (&e)e of b there exist a slow scale net (re)^ with > 1 and Sq g]0, 1] such that 

Wbeh^ix^) <\og{re) (0<e<eo). m) 

If fi is as in Theorem \3.1\ and in addition /i, /2 G tJ^oo^^g i))' '^^'^ -9 ^ ^11°° (Xt)' ^^^""^ unique 
Colombeau solution u to the initial-boundary value problem ([T]) belongs to Q^c^^^Xt)' 

Proof. From the hypotheses we may assume the following: 

(A) There exists Mq G N such that for all fc, / G N 

l|a.Vl.s|li.((o.l)) + l|5'/2.e|li2((o,l)) + Ii9^a,^<?e|li.(x.) = ©(e"*^") [s ^ 0). 

(B) For all fc, ^ G N there is a slow scale net (se)e, s^ > 0, such that 

l|5^e||L^((0,l)) + \\dld%\U^^x^j = 0{S,) (S 0). 
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(C) The constants D?p and in the energy estimate ([5^ satisfy for smaU e > (with the sfow 
scale net {r^)^ as in (|30fp ) 



0(||c, 



0(1) 



exp(F|) = 0(exp(||&,||L-)) = 0{re). 
Thus the basic energy inequahty p2p immediately yields (with Afo from property (A) above) 

= 0(r, = 0(e-^'^"-i) (e -> 0) (i?,^) 

for (/,fc)e{(0,0),(l,0),(0,l),(0,2)}. 

We proceed along the lines of Steps 1 and 2 in the proof of Theorem 13.11 and use the notation 
introduced there. 

Thanks to (A) and (B) we obtain that /i^^, /2,e, and 5^ satisfy the following variant of (A): 
(A') For all fc, ^ e N there is a slow scale net (se)e, > 0, such that 

l|5'/l,e|li.((0,l)) + l|5^72,e|li.((04)) + \\dld%\\l.(x^) = 0{s, ■ {s 0). 

When I = 1 and fc = 0, 1, 2, or / = 2 and fc = 0, the energy estimate now implies 



which yields a first extension of (Rtk I to the cases {I, k) = (1, 1), (1, 2), (2, 0). The same kind of 

observation can be repeated for the next iteration and checked explicitly in terms of /2,e, and 
cje to obtain the estimate (Ri,k I for the pairs (/, fc) — (3, 0), (2, 1), (2, 2) in addition. Successively, 
we obtain that 



(Rik) holds for alU S N and fc = 0, 1, 2. 



It is crucial to note that all the estimates (Ri^k) established so far hold with stronger upper 
bounds of the form "big oh of some slow scale net times £~^^°" and should be used in this form 
throughout the iterative process. It is only for extraction of regularity information a posteriori, 
where we deduce from these the (slightly weaker) uniform upper bound in the form 0(e~^^"~^). 
Keeping this in mind an inspection of Step 2 in the proof of Theorem 13.11 enables us to reach 
the following conclusion: Let Mg be as in property (A) above, then for all A; G N the asymtotic 
estimate 

is valid. Therefore u ~ [{ue)e] belongs to ^^ooj-j^^-j, i.e. is a regular Colombeau generalized 
function. ■ 
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